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Abstract. Measurements of ultra heavy nuclei at GeV/n energies in the galactic cosmic radiation address the question of 
the sources (nucleosynthetic s- and r-processes). As such, the determination of CR source abundances is a promising way to 
discriminate between existing nucleosynthesis models. For primary species (nuclei present and accelerated at sources), it is 
\ generally assumed that the relative propagated abundances, if they are close in mass, are not too dilferent from their relative 

^ . source abundances. Besides, the range of the connection factor associated to propagation has been estimated in weighted slab 

Ijr^ ' models only. Heavy CRs that are detected near Earth were accelerated from regions that are closer to us than were the light 

nuclei. Hence, the geometry of sources in the Solar neighbourhood, and as equally important, the geometry of gas in the same 
' region, must be taken into account. In this paper, a two zone diffusion model is used, and as was previously investigated for 

' radioactive species, we report here on the impact of the local interstellar medium (LISM) feature (under-dense medium over 

a scale ~ 100 pc) on primary and secondary stable nuclei propagated abundances. Going down to Fe nuclei, the connection 
between heavy and light abundances is also inspected. A general trend is found that decreases the UHCR source abundances 
relative to the HCR ones. This could have an impact on the level of r-process required to reproduce the data. 

CIh' Key words. Diffusion - ISM; Cosmic Rays, Bubble - Galaxy; Solar Neighbourhood, UHCR abundances 
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1. Introduction cleosynthesis) and exact abun dances are still debated ( iGorielvl 

ll999HCowan & Snederj2004 and CRs could provide some an- 

; Beyond the iron peak, the flux of Cosmic Ray nuclei drops ^^^^^ However, to determine the CR source abundances, the 

K> ■ ^^^^""^^ °^ magnitude. The UHCRs require pecuUar g^^g^ measured near Earth must be propagated back through 

■ environments to be nucleosynthesised and two distinct pro- ^^e Galaxy, which is not a straightforward task, 
cesses are involved: very generally, either the neutron flux in 

the medium is low so that flie neutronic capture rate is less than We shall not discuss here the important question of the time 
the /3-decay rate of nuclei, or the flux is high (capture rate much elapsed since nucleosynthesis and propa gation which is studied 
greater than /3-decay rate); these two situations are commonly through radioactive heavy elements (iThielemann et al.ll2002h 
referred to as s- (slow) and r- process ( rapid) . The situation is and which is another part of the puzzle, nor shall we use any 
obviously more complex - seelMeyeil \l994 for an illuminat- propagation network to study all the abundances jLetaw et alJ 
ing review -. One particular issue is the determination of the .1984) . We instead focus on propagation. The goal of this pa- 
relative contributions of s- and r-process nucleosynthesis re- pgr is to demonstrate the importance of flie local CR-source 
quired to explain the data, inferring the corresponding source distribution and in particular to examine the effect of the lo- 
abundances from the elemental a nd isotopic abu ndances. cal gas distribution on calculations of source abundances from 
Based on model predi ctions llMevei[ll994 as well as on data. As a matter of fact, whereas the truncation of path lengths 
some analysis of CR data jBrewsteretalHigO) . it is shown, for UHCRs has been discussed (through the weighted slab ap- 
for example, that for Z > 89, the entire Solar System abun- proach) and recogn ised as having a large impact on the prop- 
dances must be attributed to the r-proce ss. On the opposite , the agatio n of UHCRs Jciinton & WaddingtonI 19931 IWaddingtoiil 
latter does not contribute to Z < 40 jBinns et alJfl98l . In- |l996l), we illustrate here how this truncation is realised as a 
between, most of the elements exist as a mixture of r- and s- deficit in nearby sources in the two zone diffusion model. The 
contributions. Their origin (e.g., massive stars, explosive nu- low density gas in the local interstellar medium (LISM) is also 

of particular importance. It is shown that this under-dense le- 

Send ojfprint requests to: C. Combet, e-mail; combetOcp . dias . ie gion leads to results different from both the standard two zone 
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diffusion model (which is equivalent to a simple leaky box 
model) and the truncated weighted slab. 

The main results of this study are: i) around a few GeV/n, 
the relative abundance of Z > 60 nuclei to the lighter ones is 
decreased when the LISM features are taken into account - the 
strength of the effect is inversely related to the value of the 
diffusion coefficient (i.e. to the energy) and is directly related 
to the destruction cross section values (i.e. to Z) -; ii) the shape 
of the spectra at low energies is sensitive to the value of the 
diffusion coefficient, providing a way to determine the l atter as 
it bre aks the propagation parameter degeneracy ( Mauri iTet al.l 
observed in standard diffusion model (i.e., with constant 
gas and source disk distribution); iii) in the UHCR context of 
relative abundance determination (Pb/Pt and Actinides/Pt), a 
gas sub-density has a major effect for mixed species (species 
with both primary and secondary contributions) only. 

In §2, we give some simple arguments which allow us to 
understand why the local sub-density is expected to have a 
major influence on heavy CRs at GeV/n energies. In §3, the 
model taking into account the sub-density as a circular hole is 
exposed. Profiles in the LISM cavity as well as spectra for pri- 
mary and secondary-like nuclei under various source geometry 
assumptions are compared to those obtained in a standard "no- 
hole" diffusion model. The correction factors used to determine 
the source abundances are then derived and consequences for 
the data are discussed in §4. We then conclude and comment 
on further developments demanded by UHCRs propagation. 

2. Why a local sub-density may affect the 
propagation 

It has been previously recognised than very heavy nuclei have 
a peculiar propagation history compared to light nuclei. This 
is related to their large destruction cross sections which make 
them have very short path lengths. Regarding this extreme sen- 
sitivity to nuclear destruction, an UHCR breaks-up more often 
than a lighter one and thus propagates for a shorter distance. 
lKaiseretalJ ( ll972h emphasized that the strength of the source 
and its location in space and time may be a dominant local con- 
tribution or merely adds to the general background depending 
on whether this source is close or far. These authors use a very 
simple model in which propagation starts from a single source 
in space and time. They conclude that A > 8 1 elements cannot 
be explained by any source older than 10^ yr In the present 
paper, the steady-state approximation is however assumed but 
the spatial discreteness hypothesis is relaxed to some extent. 
This is a first step away from a leaky box model towards a 
more realistic description. The framework is a two zone diffu- 
sion model (thin disk and h alo) which has been used by many 
authors for light nuclei, e.g. iM aurin et al]ll200lll . In particular, 
such a model has distinct features compared to leaky boxes. 
The interstellar matter density is now located in a thin disk 
with CRs spending most of their time in the diffusive halo. The 
characteristic times of the various processes i n competition dur- 
ing pr opagation (see below) can be extracted jTaillet & MaurinI 
As it is shown below, the fact that the typical distances 
travelled by heavy or UHCRs before their destruction could be 



as little as a few hundreds of pc naturally means that LISM 
properties are an important ingredient for propagation. 

2.1. Propagation characteristic times vs A 

During the propagation of nuclei through the Galaxy, sev- 
eral processes affect their journeys. All of these have specific 
time/length scales which characterise the importance of the 
considered process relative to the others. The nuclei can un- 
dergo energy losses of different forms, interact with atoms of 
the ISM and be destroyed in a spallation process or escape from 
the boundaries of the diffusive volume (i.e. the galactic halo) 
because of the combined influences of diffusion and convec- 
tion. 

In a leaky box model, the density of gas in the box is con- 
stant so that the link between characteristic lengths and times 
is simple and the comparison to the escape length straightfor- 
ward. In reality, spallation and energy losses occur when a nu- 
cleus crosses the galactic disk. Thus, the characteristic times (or 
lengths) for those pro cesses are directly link ed with the num- 
ber of disk-crossings dTaillet & Maurinl2003h which should be 
compared to the diffusive and convective escape times. 

2.1 .1 . Spallation dominates over energy losses 

Whether in a leaky box model or a more sophisticated model 
as used here, spallations and losses both occur when gas is tra- 
versed. We assume its density to be nisM = 1 cm""*. The spalla- 
tion rate is given by Fjp = nisMO"v, where v is the velocity of the 
nucleus and cr is th e reaction cross section. For this study, the 
iLetaw et al ] (ll983l) cross sections are accurate enough (Ek/„ is 
the kinetic energy per nucleon): 

o- = cr^[l - 0.62 exp(-£^/„/200) sin(10.9£^°„2'^)] , 

o-oo(Ek/n > 2GeV/n) = 45A'"[1 + 0.016 sin(5.3 - 2.63 log(A))] 

The loss rate is given by Hoss = {dE I dt\g^^ I Ek, 
where Ek is the kinetic energy of the nucleus, and the 
Coulomb and ionisation energy losses are taken into accoun t 
jMannheim & Schlickeise r' 1994: ' Strong & Moskalenkol 19981) 
assuming the ionised fraction of the ISM to be 0.033. In both 
cases, dE/dt is proportional to Z^, with Z the charge of the 
nucleus. 

At a given Ek/n, Fsp scales with the atomic mass A as A"'^ 
(log A does not induce strong variations). Assuming A - 2Z, 
which is not the case for heavy nuclei (A > 2Z), ionisation 
and Coulomb loss rates also increase with A but following 
Tion/coui A. A comparison of these two rates yields at 1 GeV/n 

F,p~ 10x(A/10)-°-^Fion, 
F,p~60x(A/10)-°-^F,oui. 

This means that the effect of spallation is always dominant over 
energy losses (at higher energy the effect is even stronger). The 
latter are thus discarded in the rest of this paper as it mainly 
focuses on qualitative effects. Note also that the cross sections 
have not been corrected for the effects of decay which can be 
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Fig. 1. Characteristic times of diffusive/convective escape and 
spallation as a function of the atomic mass. Typical values 
Kq - 0.0112 kpc-Myr"' and Lhaio - 4 kpc were considered. 
The times have been calculated at kinetic energies per nucleon 
of 1 GeV/n and 10 GeV/n in the two zone diffusion model. 



significant for the propagation of heavy nuclei. These effects, as 
well as energy losses and reacceleration, will be properly taken 
into account in a separate paper where y6-decay and electronic 
capture decay will be re-examined in the context of the present 
diffusion models. 

2.1 .2. Time scales of diffusion, convection and 
spallation 

Once established that spallations are dominant relative to the 
energy losses, we compare the former to the convection and 
diffusion processes taking into account the number of disk- 
crossings in the two zone diffusion model. For further justi- 
fication and details the reader is referred to iTaillet & MaurirJ 
1I2OO3I) . 

It is well known that rjiff = L^/K{E) is the characteristic 
escape time needed by a nucleus to leave a diffusive volume 
of size L for a diffusion coefficient K{E). The presence of a 
galactic wind (velocity Vc), assumed to be constant and per- 
pendicular to the galactic plane, induces a general convective 
motion throughout the Galaxy. Considering that the diffusive 
process is still present, the boundary to be taken into account is 
not the halo size but the distance for which the diffusion can- 
not compete anymore with the convection and brings the nu- 
cleus towards the Earth. With this consideration, one can de- 
fine the typical time scale for the convection escape process as 
Twind = 2K(E)/Vl 

Spallation only occurs within the disk of total thickness 2h. 
The typical length associated with the spallation process in a 
diffusive propagation is r^j, - K{E)/{hrsp). The time scale, as- 
suming a purely diffusive transport is then r^p - r^^/K{E) = 
K(E)/(hT,p)\ 

The characteristic times for diffusion, convection and spal- 
lation are plotted on Fig. ^ as a function of the atomic mass 
of the nucleus A at two different kinetic energies per nu- 



cleon (1 GeV/n and 10 GeV/n). As we assume K{E) — 
KoPCR/lGeWy (where is the rigidity), the convection and 
diffusion times are constant with A. At 1 GeV/n, it appears 
that spallation is almost always dominant, except for the very 
light nuclei where both convection and diffusion compete. At 
10 GeV/n, spallation is less dominant but remains the major 
process for the heavier nuclei. A naive (and false) interpreta- 
tion of these figures is that a heavy nucleus never escapes from 
the Galaxy. Actually, many paths lead from the source to the 
halo boundary and these numbers only indicate the typical time 
a CR reaching our position can have travelled before being de- 
stroyed. For convection and diffusion, Tdifr and r„ind define an 
exponential cu toff whereas the cutoff is smoother for spalla- 
tion (Taillet & MaurirJ2003l) . 

Tabled contains the typical spallation distances for differ- 
ent atomic masses (10, 50, 100 and 200) and a set of three dif- 
fusion coefficients (two extreme values and one median) which 
fit the B/C ratio when combined with t he appropriate choice fo r 
the other propagation parameters OMaurin et alJl200ll l2002h . 
For A=100 and A=200 and for the median value of the dif- 
fusion coefficient, the spallation length scales have the same 
order of magnitude as the size of the observed local gas sub- 
density (see Sect.O. This sub-density precludes the creation of 
secondary species in the local bubble and the average density 
crossed by a primary or a secondary nucleus during its journey 
is smaller than when this sub-density is not taken into account. 
It is expected that heavier nuclei will be more sensitive than 
lighter nuclei to this feature. 



A 


Group 




Ko (kpc^ Myr"') 








0.0016 


0.0112 


0.0765 


10 


LiBeB/C 


210 pc 


1.25 kpc 


6.38 kpc 


50 


Sub-Fe/Fe 


70 pc 


390 kpc 


2.01 kpc 


100 


Z=44-48 


40 pc 


250 pc 


1.27 kpc 


200 


Actinides 


30 pc 


150 pc 


780 pc 



Table 1. Typical length scales of the spallation process for dif- 
ferent nuclei and diffusion coefficients at 1 GeV/n. 



The smaller the spallation length scale, the greater the effect 
of the sub-density is expected to be and this suggests that the lo- 
cal bubble must be considered to treat the heavier species. This 
reasoning is valid only if the cosmic rays undergo the same dif- 
fusive process in the LISM than in the rest of the Galaxy. There 
is at least two indications supporting this assumption: 

- The local bubble and the galactic halo share some common 
properties. Hence, it seems reasonnable that the values of 
the diffusion coefficient in both regions are not drastically 
different. As we use only one coefficient for the disk and 
the halo, the latter can be seen as an effective value that 
also applies for the local bubble. 

- The second clue comes from the analysis of secondary ra- 
dioactive species which can travel for a few hundreds of 
pc only before decaying. tPtuskin & SoutouH ( Il998l) were 
able to derive the corresponding 'local' diffusion coefficient 
modelling the LISM as three shells of gas with various den- 
sities. They found a larger value at low energy than the 
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standard one found from B/C analysis. This again points 
towards a diffusive transport of CRs in the local bubble. 

3. Two zone model with a hole 

Firstly, the local interstellar medium (LISM) is known to be a 
highly asymmetric low density region, extending between 65 
to 250 pc - recent studi es have for the first tim e conducted to- 
mography of the LISM jLallement et al.lE003l) -. As a conse- 
quence, the spallation rate is lessened so that the primary fluxes 
are expected to be enhanced. 

S econdly, the history o f the origin of the local bubble 
(see iBreitschwerdt & Coxl ( l2004l) for a review) is an im- 
print of t he explosive steflar a c tivity in the Solar neigh- 
bourhood dMaiz-ApellanizlEoOlt iBerghofer & Breitschwerdll 
l2002t iBenitez et alJ l2002h . The spatio-temporal position of 
these sources is certainly of importance. However, as the 
steady state is assumed in this paper, o nly the spatial in- 
fluenc e can b e inspected. For example, iLezniak & WebbeJ 
(I1979I) : Iwebber. (.1993i) have shown that an absence of nearby 
sources (in a diffusion model) can truncate the path lengths 
in weighted slab models. This truncation has been recognised 
to b e of great importance in the propagatiori of h eavv nu- 
clei dChnton & Waddingtonll993tlWaddingtonll996l) . This ef- 
fect is realized here as a circular hole of a few hundreds of pc 
surrounding the Solar area reflecting the fact that no very recent 
sources have exploded in the Solar neighbourhood. Note that 
this truncation could also be an effect of the matter traversed 
during the acceleration and the escape from a source region. 

3.1. The model 

The effects of a local sub-density on radioactive species 
have been studied in the context of a diffusion model by 
iDonato et alJ (l2002l) . As it is difficult to derive an analytical 
solution of the diffusion equation using a realistic distribution 
for the gas, it is assumed that the LISM is a circular cavity of 
radius a (see below). It appeared that the radioactive flux re- 
ceived on Earth was lessened by ~35% when the decay length 
of the radioactive species' was twice the size of the sub-density 
compared t o the case wh e n a su b-densitv was not taken into 
account (cf. IPonato et alJ ll2002i) . Fig. 4). Furthermore, in or- 
der to match measured radioactive abundance ratios ( "'Be /''Be, 
^''Al/^^Al and ^''Cl/Cl), the size of the local sub-density was 
constrained to lie within 60 and 100 pc - values consistent with 
direct observations. This increases confidence in the model. In 
the latter, the Galaxy is considered as a thin disk where the 
sources and gas are located (Fig. |2j- T o be consistent with 
the results obtained for light stable nuclei ('Ma urin et all200ll) . 
for radio active nuclei (iDonato et al.i .2002) and also for an- 
tiprotons dPonato et al.l200ll) . the same geometry is used. The 
width of the disk, 2h, is 200 pc and the Galactic radius R is 
20 kpc. The exact value for R is not important so that one can 
set the center of the cylindrical geometry at the Earth posi- 

' In a diffusion model, the typical length for decay is defined as 
'rad = ^Kjtq where K is the diffusion coefficient, 7 the Lorentz factor 
and To the lifetime of the nucleus. 




L 



Fig. 2. Schematic view of the cylindrical geometry for the 
model used in this paper: the hole region of radius a is located 
in a thin disk surrounded by a diffusive halo of size L (a convec- 
tive wind may be present). The exact radius of the Galaxy 
radius R is not important (see text). 



tion (see iDonato et alJ ( l2002i) ). which simplifies the calcula- 
tions. The diffusive halo, empty of gas and sources, extends to 
a height L on each side of the disk which can vary from 1 to 
10 kpc. The local sub-density is then a hole of radius a in this 
disk, whereas the source sub-density is a hole of radius b. From 
the disk emerges a galactic wind with a constant velocity Vc 
that adds a convective component to the diffusive propagation 
of the CRs. 

These assumptions (sources and gas in a thin disk, circular 
hole, halo empty of gas, same diffusion coefficients in the disk 
and in the halo) are probably too strong. However, they have 
lead to a successful description of many species in a consis- 
tant way, an issue we think to be of great importance. We also 
believe that the approximations made (some that would need 
a heavy numerical treatment to be relaxed) would have a mi- 
nor effect compared to the one induced by the LISM geometry. 
Despite many weaknesses, the present diffusion model has to 
be thought as a first step towards a more realistic description. 

Energy losses are not taken into account meaning that in 
this cylindrical geometry the diffusion equation reads 

- K^N +Vc^ + 2hT6{z)N = q(r, z)6{z) (1) 
oz 

where - N{r,z) is the differential density in energy. The 
Dirac distribution 6{z) expresses the fact that the sources and 
the gas are only located in the thin disk. In this work, the re- 
sults are always normalized to the source spectrum thus the 
only quantity one must specify is q{r), the radial source dis- 
tribution (cf. lA.H . It is appropriate to use a decomposition in 
Bessel space to solve this equation and the complete derivation 
is given in d etail in ApperidixlXI( the derivation closely follows 
that given in lDonato et alJ ( l2002i) ). 

The height L of the halo is a free parameter as are the galac- 
tic wind velocity Vc, the spectral index 6, and normalization Kq 
of the diffusion coefficient. Those parameters must be tuned to 
fit the data. In the following work and when mentioned, the 
max, median and min set of parameters can be understood as 
detailed in Table |2 These three sets are the two extreme and 
median sets (with regard to the value of the diffusion coeffi- 
cient normalization) shown to be compatible with B/C analy- 
sis. Eventually, throughout the paper, the gas sub-density is set 
to fl = 100 pc whereas the sub-density in sources can be varied. 
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set 


6 


Ko (kpc- Myr') 


L (kpc) 


V, (km s-') 


max 


0.46 


0.0765 


15 


5 


med 


0.7 


0.0112 


4 


12 


min 


0.85 


0.0016 


1 


13.5 



Table 2. Two extreme and one m edian sets of param eters 
shown compatible with B/C analysis jMaurin et al.l200ll) . 



In the remainder of this section, we focus on the differences 
between a model with a hole in the surrounding gas (and/or 
sources) and the standard diffusion model with no hole. Note 
that as the propagation parameters used for the standard diffu- 
sion model (with no hole) are fitted to B/C, the results obtained 
in this latter model and those obtained in a LB are similar (see 
Fig. |8}. It is useful to define the "enhancement factor" which 
is defined as the ratio of a given hole configuration to the stan- 
dard diffusion model (no hole, hereafter SDM) for a given set 
of propagation parameters. Such a ratio ensures that the source 
spectral shape and normalization are factored out for primaries 
and the production cross section is for secondaries. The discus- 
sion of the absolute effects on abundances - including the stan- 
dard result from the leaky box model - is contained in Sect.|3 

3.2. General behaviour 

To separate the effects of the absence of gas and sources in the 
hole, we study three different configurations: /) in the first case, 
only the gas is absent from the "hole region" but the sources are 
still present, /;) the second case is the opposite configuration 
and the hole is a sub-density of sources but not of gas and ;;;) 
finally, we consider the situation where there is both a hole in 
gas and a hole in sources (possibly having different sizes). 

3.2.1. Profiles 

The enhancement factors for primaries, secondaries and for the 
ratios Secondary /Primary are plotted on Fig.|3as a function of 
the distance to the center r (r = corresponds to Earth lo- 
cation). The three hole configurations described above have 
been studied. For these profiles, the size of the hole is 0.1 
kpc, and the median set of parameters has been used, namely: 
Ko = 0.0112 kpc^ Myr-', L = 4 kpc, Vc = 12 km s"' and 
6 - 0.7. The iron group (left panel) and the actinide group 
(right panel) are both considered at 1 GeV/n. Firstly, it appears 
that when the distance from the hole is large enough (typically 
3a), the enhancement factor tends to unity (i.e. there is no dif- 
ference from the no-hole case) emphasing that, in a first ap- 
proximation, the hole has only a local effect. An important con- 
sequence is that nuclei are not sensitive to other density features 
of the ISM, which justifies the crude model used here. 

When a hole in gas only is considered, more primaries (up- 
per panels) are expected (i.e., enhancement factor > 1) com- 
pared to the SDM as nuclei entering the hole region do not un- 
dergo spallation as they would if gas were present. In that case, 
at the Earth location, there is a 60% increase in primaries of the 
iron group and a 275% increase in the actinides. In the oppo- 
site configuration (sources absent, gas present), no primary nu- 
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Fig. 3. Left panels: Fe-group. Right panels: Actinide group. 
Primaries (upper), secondaries (middle) and sec/prim (lower) 
enhancement factors are plotted as a function of the distance 
to the center of the sub-density (the fluxes are normalized rel- 
atively to the SDM fluxes). In both models, the hole in gas 
and the hole in sources are equals a - b = Q.l kpc and 
£'jt/n=l GeV/n. The set of astrophysical parameters is the me- 
dian set: Ko = 0.01 12 kpc^ Myr"', L = 4 kpc, V,. = 12 km s"' 
and 6 = 0.7. 



clei are produced within the sub-density and those which have 
propagated to the hole undergo spallation. In that case, the pri- 
mary density is lessened compared to the SDM (~ 20% for 
Fe-group nuclei and ~ 40% for actinides). 

The third case, where both gas and sources are absent from 
the hole region is found to give intermediate results. It appears 
that in this geometry the absence of gas has a larger impact 
that the absence of sources: for a hole in gas and sources, the 
enhancement factor is greater than unity. There is a 20% ef- 
fect for nuclei of the iron group and a 75% effect for actinides. 
However, the result is quite sensitive to the particular geometry 
assumed. To obtain a very rough estimate of this sensitivity, the 
diffusion equation was solved for a spherical geometry where 
sources and gas holes were taken as shells. In this latter config- 
uration (which is unreaHstic, but which may be considered as 
the extreme opposite geometry to that of the thin disk case) the 
enhancement and relative importance of a hole in gas or a hole 
in sources are smaller and the gas-hole enhancement is almost 
cancelled by a hole in sources. 

When looking at the secondary enhancement factor (Fig.|3l 
middle panels), a sub-density in sources only has almost no 
consequences (~ 10% effect) on the secondary flux (enhance- 
ment factor ~ 1). This indicates that most of the secondaries 
found in the solar neighbourhood originate from primaries that 
were not produced locally. When a gas hole is considered (close 
sources are present), the secondary enhancement factor is max- 
imal (~ 80% for Fe-group nuclei and ~ 450 % for actinides). 
In this case, the secondaries present at r = cannot have 
been produced locally (as there is no gas) and have all been 
propagated from further regions. Once they reach the gas sub- 
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density, they cannot be destroyed, which explains the enhance- 
ment. The enhancement in the case of a hole in gas and sources 
is explained in the same manner. However, this enhancement 
is slightly lower (~ 60% for Fe-group nuclei and ~ 350% for 
actinides) than the case in which only gas is excluded: when 
sources are present in the hole, primaries from the hole can 
propagate in the disk and produce secondaries that may diffuse 
back towards the Earth. 



3.2.2. Spectra 

Enhancement factor spectra for primaries, secondaries (Fe- 
group nuclei in the left panel and actinides in the right panel) 
and for the ratio Sec/Prim are plotted in Fig. 0] between 
500 MeV/n and 100 GeV/n at r = 0. At high energy the be- 
haviour of any hole configuration tends, as expected, to the no- 
hole case at high energies as spallation becomes negligible. In 
the case of a sub-density in gas, the enhancement factor for the 
median propagation configuration is quite large (almost a fac- 
tor of 2), even for the Fe flux. Comparatively, the sub-Fe/Fe 
ratio is enhanced by a mere 20%. 



Fe-group 



Actinide group 




Fig. 4. Same as Fig.|3but for spectra at r = 0. 



Actually, these enhancements, whether for primaries or sec- 
ondaries, are not very sensitive to the exact value of L and Vc 
(because the enhancements are normalized to the SDM). They 
mostly depend on the value of the diffusion coefficient. Hence, 
it implies, for example, that a very precise measurement of the 
iron abundance spectrum would give some hint as to the value 
of the diffusion coefficient, by evaluating the deviation from a 
standard diffusion model. The larger the deviation at low en- 
ergy, the smaller the diffusion coefficient required to fit the de- 
viation. 



3.3. Other dependences 
3.3.1 . Enhancement factor vs 

It is useful now to display the enhancement factor as a function 
of fl/ (fl being the size of the hole and r the typical spalla- 
tion distance, see Sect. 12. 1.21 1. As several co mbinations of (Kg, 
6) for a given choice of A yield the same jTaillet & MaurinI 
I2OO3I) . the parameter dependence is more economically stud- 
ied. Indeed, to quickly obtain the enhancement due to a hole 
in gas - or any combination of holes, compared to the SDM - 
it is sufficient to take the desired Kq and 6, choose a specific 
A and evaluate for a given energy. The enhancement fac- 
tor is then directly inferred from Fig. |5l Two cases are dis- 
played to emphasize the effect of the Galactic wind: the first 
with Vc - km s ' and the second with Vc - 12 km s"'. This 
is the typical range of the possible values for the wind. It is not 
a dominant effect. 



Primaries 



I • I r 
gas and source hole 
gas hole 
source hole 



Thiols lines: Vc-0 km/s 
Thin lines: Vc-12 km/s 




0.5 



Fig. 5. Primary and secondary enhancement factors as a func- 
tion of the ratio a /rsp for the median set of parameters. The 
size of the holes are a = b - O.l kpc, and the different values 
of cover a wide range of energies, Kq and destruction cross 
section values. 



The enhancement factor for several nuclei could be ob- 
tained from the previous plot, but for simplicity, they have been 
gathered for light to heavy nuclei for several energies and prop- 
agation sets in Table |3] Despite the fact that the three sets of 
propagation parameters have been shown to be compatible with 
B/C analysis, the minimal set gives some unrealistic enhance- 
ment factors. This allows the exclusion of an unrealistic config- 
uration fo r which the diffusion slope required to fit the B/C data 
was 0.85 jMaurin et al.ll2002l) . far away from any theoretical 
expectations. Actually, in previous papers attempting to obtain 
some conservative value for the pr opagation parameters from 
the sole B/C and sub-Fe/Fe ratios jMaurin et alJl200lLl2003) . 
the hole configuration has not been taken into account. In or- 
der to be able to firmly exclude some Kq and 6 values, a new 
detailed analysis is required, which goes beyond the scope of 
this paper. In conclusion, it seems that very small values of Kq 
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have to be excluded because of the too large enhancement they 
yield when a hole in gas is included. 





Prim 


Sec 




min/med/max 


min/med/max 


LiBeB-CNO 






E=l GeV/n 


2.08/1.17/1.03 


2.29/1.13/1.02 


E=5 GeV/n 


1.36/1.07/1.02 


1.33/1.05/1.01 


E=10GeV/n 


1.21/1.04/1.02 


1.19/1.03/1.01 


Fe-group 






E=l GeV/n 


6.24/1.58/1.11 


20.9/1.75/1.11 


E=5 GeV/n 


2.34/1.23/1.06 


3.38/1.26/1.06 


E=10GeV/n 


1.71/1.15/1.04 


2.03/1.16/1.04 


Z=44-48 






E=l GeV/n 


8.67/1.80/1.14 


43.0/2.14/1.15 


E=5 GeV/n 


2.95/1.31/1.08 


5.06/1.36/1.08 


E=10GeV/n 


1.99/1.19/1.06 


2.59/1.22/1.06 


Actinides 






E=l GeV/n 


17.6/2.77/1.26 


200/4.45/1.31 


E=5 GeV/n 


5.73/1.61/1.15 


18.2/1.82/1.16 


E=10GeV/n 


3.28/1.37/1.11 


6.22/1.46/1.12 



Table 3. Enhancement factors at the Earth's location due to 
the presence of a local sub-density on four groups of nuclei. 
Different energies (from 500 MeV/n to 10 GeV/n) and three 
sets of parameters have been considered. 



' 1 ' 1 


1 1 1 1 1 




— Primaries, Fe-group 




— Primaries, Actinide group 




— Secondaries, Actinide group 




— Secondaries, Fe-group 






1 , 1 




1 7^^~~i"""""' 



b (kpc) 



Fig. 6. Enhancement factor for primaries and secondaries of 
the Fe and actinide groups as a function of the size of the hole 
in sources, b. A constant gas hole radius a = 100 pc was con- 
sidered and the propagation computed using the median set of 
parameters. 



nearby sources c annot be easily dis connected to the question 
of recent sources jTaillet et alj2004 . 



3.3.2. Impact of the hole in sources b 

A few decades ago, some authors introduced the possibility of 
a truncation of th e path lengths us ed in the weighted slab for- 
mahsm. [Lezniak & Webbe J (ll979h then showed that a hole in 
sources in a diffusion model mimics such path-length trunca- 
tion. This truncation was first introduced to fit best the sub- 
Fe/Fe ratio, though afterwards, new results on cross section 
production reduce the need for it. However, truncation was 
implemented for UHCR pr opagati on and it has a sizea ble ef- 
fect (Clinton & Waddingtonl ll993l: IWaddingtonl ll996l) . espe- 
cially for secondary nucleus production. The truncation uncer- 
tainty is here equivalent to the uncertainty in b (see also Fig.|8] 
in Sec. O- 

In this section, the hole in sources varies from b - to 
b — 500 pc. In Fig.|6l the enhancement factors of primaries and 
secondaries of the iron group and actinide group are plotted as 
function of the size of the hole in sources. The energy of the 
nuclei is 1 GeV/n and the size of the gas hole is once again 
fixed at fl = 100 pc (because this feature is certainly present). 
Depending on the hole size b, one can have an enhancement 
factor that is greater or lower than unity. Actually, this situation 
occurs only because the hole in sources is in competition with 
the hole in gas. 

However, we do not wish to put too much emphasis on this 
effect. It must be considered as an illustration of the impact 
on relaxing the condition of continuous distribution of sources. 
The limits of our model are reached and to go further, one must 
relax the steady state approximation. In fact, the question of 



3.4. Preliminary conclusions 

Setting holes in a diffusion model affects the abundance spec- 
tra, especially at low energies where the importance of spalla- 
tion is larger A hole in gas increases the number of primaries, 
even more the number of secondaries. Conversely, a hole in 
sources decreases the number of primaries (the secondaries be- 
ing quite insensitive to it). So there is a balance between the 
two configurations. Unless a large hole in sources is chosen 
{b > 300 pc), because the hole in gas (size a ~ 100 pc) 
is firmly established - from direct observations, and also in- 
directl y because it better fits light radioactive CRs measure- 
ments ( iDonato et al .l2002h -■ the final effect is an enhancement 
of the cosmic ray fluxes. 

The dependence of this enhancement to the halo size L 
and the convective wind Vc is minor The main dependence 
is through the diffusion coefficient, K{E). The previous sec- 
tions helped us to establish that very small diffusion coefficients 
that were f ound to best fit the B/C ratio in a standard diffu- 
sion model (iMaurin et alJ2002l) must be discarded as they pro- 
duce a too large enhancement of fluxes when the hole in gas is 
taken into account. It is reassuring since these small values for 
Kq corresponded to quite large diffusion slopes, unsupported 
by theoretical considerations (too far from a Kolmogorov or 
a Kraichnan turbulence spectrum). However, a new study of 
B/C must be conducted taking into account the hole configu- 
ration to provide a quantitative result. Going one step further, 
one sees that the low energy secondary fluxes can be useful to 
determine Kq. The secondaries better suit this estimation than 
primaries as, unlike the latter, they are not very sensitive to a 
hole in sources. This could provide another way, apart from 
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using radioactive nuclei dPtuskin & Soutoullll998l) . to extract 
the diffusion coefficient without too much ambiguity. Note also 
that the heavier the nucleus, the more local it is, so that look- 
ing at different nuclei gives different sampling regions where 
the diffusion coefficient is averaged. Figure Q in the next sec- 
tion will provide a better understanding of how the effect of 
the hole sometimes disconnects from the dependence on L and 
so in certain configuration why the flux at low energy only de- 
pends on Kq (if the latter is not too large). 

One can now come to the differences in derived source 
abundances. Cosmic ray fluxes are propagated to source by us- 
ing a leaky box (or equivalently a standard diffusion model) 
and hole models applied to real data. 

4. Consequences for CRs observations 

UHCR abundance spectra have been obtained via space- 
craft me asurements s i nce t he 1970s, most notably by 
ARIEL 6 frowler et al.' 'l987|), heao-3 dBinns et all Il989l) . 
TREK fwestph al et al. 1998) and the Ultra-Heavy Cosmic-Ray 
Experiment (or uhcre; Donnelly et al., in preparation). 

The data are scarce, especially in the actinide region, and 
only elemental abundances have been obtained. However, sev- 
eral important conclusions have already been drawn. They are 
related i) to nucleosynthesic aspects (see the introduction), ii) 
to the possible site of acceleration and iii) to the mechanisms 
leading to elemental segregation during acceleration. 

The recent results of uhcre provide unprecedented statistics 
and will allow to give firmer conclusions. However, as empha- 
sized throughout this paper, some aspects of the propagation 
still need to be clarified in order to fully interpret these data. 
In this last section, some questions and results about UHCRs 
are first reiterated, before inspecting the propagation effects on 
these nuclei in different models. In particular, those species 
suffering from major propagation uncertainties (regardless of 
cross-sections, which is another issue) are sorted. 

4.1. Introduction: UHCR data, their interpretation and 
general issues 

The ultra-heavy (Z > 70) elemental abundance ratios most per- 
tinent in determining cosmic ray origin are Pb-group/Pt-group, 
92U/9oTh and Actinides/Pt-group^. 

The first, Pb/Pt, provides clues as to biases in the CR-source 
abundances. Elements in the Pt-group are mainly intermediate- 
FIP refractory and r-process, while in the Pb-group low-FlP, 
volatile and primarily s-process elements predominate. Thus 
an anomalous GCR Pb/Pt ratio (relative to Solar values) would 
show up any nucleosynthetic (s- or r-process) or atomic (FIP 
or volatility) bias in the source abundances. This ratio is more 
enlightening than (for example) the Pt/Fe or Pb/Fe ones, as the 
ratios of nuclei similar in mass, are supposed to be relatively 
unaffected by interstellar propagation. This point is discussed 
in the next section. In common with other space-based mea- 
surements, the UHCRE results (Donnelly et al., in preparation) 

2 Pt-group= (73.5 < Z < 80.5), Pb-group= (80.5 < Z < 83.5) and 
Actinides are Z > 88. 



demonstrate that the Pb/Pt abundance ratio is decidedly low in 
the OCRs (0.24 + 0.03) compared to the best estirnates f rom so- 
lar and meteoritic material (1 .03 + 0.12: lLoddersl ( l2003l) V Even 
assuming a very severe propagation effect on this ratio (x2.6), 
the GCR value is a mere 0.63 + 0.09. This could indicate a 
volatilit y-based accel eration bias as Pb elemen ts are relatively 
volatile jElhson et alAl997i:.Meyer et alJl99^ . 

The relative abundance of chronometric pairs such as 90 Th 
and 92U can provide an estimate of the time elapsed since 
their nucleosynthesis. Most models of actinide decay indicate 
that 92U/9oTh drops to unity about 1 Gyr after nucleosynthe- 
sis. However, only 44 actinides have been detected so far. 35 
of these were detected by the uhcre and this experiment pro- 
vides the only estimate yet of the 92U/9oTh in the CRs. The 
Icr upper limit is only ~ 1, implying that significant 92U de- 
cay has occurred and that the time elapsed since nucleosyn- 
thesis is relatively large. The transuranics can also be used as 
excellent cosmic ray 'clocks' since the relative abundances of 
q^N p, 94PU and 9fiCm fall dra stically 107 yr after nucleosynthe- 
sis jBlake & Schrammll974l) . Again, the best data available are 
from the uhcre, which detected one geCm event. The longest- 
lived isotope of Cm has a half-life of only 16 Myr. This fact, 
combined with the 92U/9oTh results from the same experiment 
suggests that the CR-source material is an admixture of old and 
freshly-nucleosynthesised matter, such as that found in super- 
bubbles. 

Finally, anomalies in the Actinides/Pt ratio could indicate 
an unusual, possibly freshly-nucleosynthesised component in 
the cosmic ray source matter Results from all of the space- 
based measurements indicate a high value relative to solar sys- 
tem material, though the uncertainties on the latter are large. 
The uhcre's Actinide/Pt ratio-' (0.028;^qqqj) in broad agree- 
ment with other observations, is higher than in the present in- 
terstellar medium (0.014 + 0.002) and similar to that of the 
protosolar medium (~ 0.023) and the interior of superbub- 
bles (0.029 + 0.005). These observational values are unadjusted 
for propagation and so are therefore best considered as lower- 
limits. 

There are large uncertainties on the effects of propagation 
on these ratios. For example, estimates of the effects of propa- 
gation on Pb/Pt vary from factors as low as xl.3 to as high as 
x2.6. Accurate estimates of propagation effects are therefore 
crucial to interpret the data. 

4.2. UHCR abundances and the correction factor for 
the propagation 

As emphasized above, in the context of UHCRs, one issue is 
the determination of relative source abundances of, for exam- 
ple, Actinides/Pt and Pb/Pt. To this end, in this subsection, it 
is verified that the choice of a more refined propagation model 
leaves the abundance ratio of close heavy nuclei (if considered 
as being pure primaries) almost unchanged. Figure shows 
that for any hole configuration and any propagation parame- 
ters, the maximum effect on the Actinides/Pt ratio (when both 
considered as pure primaries) is at most 20% (as obtained in the 

3 i.e. (Z > 88)/(75 < Z < 79). 
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source hole case). The effect is even smaller for the Pb/Pt ra- 
tio (also considered as pure primaries). However, Pt is a mixed 
species and is responsible for the uncertainties on these two 
ratios. Considering for example the ratio Pt^'^'^/Pb''"™, the dif- 
ference in Pt production from Pb in the various hole model can 
be as large as a factor of 2. Hence, this hole must be taken into 
account using a complete propagation network to derive abun- 
dances in a model that is thought to be more realistic than the 
LB model. 

Figure is also useful in understanding the propagation 
regimes of heavy nuclei. For very small diffusion coefficients 
(equivalent to low energies), the spallation length is much 
smaller than the halo size. In that case, the ratio of two close 
nuclei is independent of L but sensitive to any hole. Conversely, 
for large diffusion coefficients (equivalent to high energies), the 
spallation length is much larger than the hole size, thus, the ra- 
tio is only sensitive to L. Depending on the nucleus and en- 
ergy under consideration, one can be in an intermediate regime 
where both the hole and halo influence the ratio. 



0.95 




0.85 



gas and source hole 

gas hole 

source hole 

no hole 



0.001 



0.01 



0.1 



Ko (kpc /Myr) 



Fig. 7. {Actinides/Pt)pi-opag /{Actinides/Pt)source ratio (all 
species considered as pure primaries). For large diffusion 
coefficients, the effect of the hole (any configuration) is negli- 
gible and only depends on the halo size L: the curves with the 
same L tend to the same values - upper (resp. lower) group of 
curves: L - 2 kpc (resp. L = 10 kpc). On the other hand, for 
low diffusion coefficients, the ratio is completely insensitive to 
L. See text for details. 



4.3. Abundance normalization bias vs A 

The ARIEL 6 experiment dFowler et alJ measured the 

abundances of Fe-group elements so that they provided abun- 
dances normalized to Fe. One peculiar feature is the overabun- 
dance of the 44 < Z < 48, 62 < Z < 69 and 70 < Z < 74 
groups which are presumed to be predominantly secondary in 
origin. Actinides w ere also found to be overabundant in this ex- 
periment. Later, in iBinns et al.l ( Il989l) . CR abundances (from 
HEAO 3) were propagated back using a LB model. This study 



rises the possibility that a bias occurs during the process be- 
cause the propagation correction factor in a LB may be incor- 
rect. 
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Fig. 8. Propagation correction factor for a given A relative 
to the nucleus A = 50. From the curves, one directly gets 
^Modei 1/^Modei 2 ^ Corr.Factor^""*"' ^/CoiT.Factor^"''"' ' which 
is the relative CR abundances obtained using different propa- 
gation models. The standard DM and the LB model (2 dash- 
dotted lines) are equivalent. The hole source DM and the TWS 
(2 dotted lines) which predict a smaller primary density than 
the LB are also equivalent. Left panel is for primaries, right 
panel for secondaries (it is assumed for simplicity that these 
secondaries come from the parent Aprin, - Asec + 6). See text 
for discussion. 



Figure|Slgives the correction factor according to the propa- 
gation scheme used. As before, for primaries, it corresponds to 
the flux divided by the source term whereas for secondaries, it 
corresponds to the flux divided by the production cross section. 
All the results are normalized to the iron group (A ~ 50). For 
the diffusion models (without hole and with different hole con- 
figurations), the median set of parameters is used {Ko = 0.01 12 
kpc^ Myr '). We also plotted the results of a LB model and 
for a truncated weighted slab (TWS). The densities in the TWS 
are calculated using the escape mean free path and pat h length 
distribution (PLD) from lClinton & Waddingtonllll993h - detail 
can be fo und in appendix 151 The truncation was taken to be 
1 g cm'"^: Iwaddingtoiil ( 1 19961) showed that such a truncation 
combined with a pure r-source better fitted the UHCR abun- 
dances. 

There is a general trend showing that the propagation in a 
LB (or as well as in a standard diffusion model) predicts larger 
fragmentation of nuclei than a diffusion model with a hole in 
gas, especially for high A. For a larger diffusion coefficient, this 
effect would be weaker Note that the effect is stronger for sec- 
ondaries (right panel) than for primaries (left panel), as already 
discussed. Furthermore, it has to be noticed that both the TWS 
model and the source hole diffusion model give lower primary 
densities that a LB or a standard DM model. This is pointing 



10 



Combet et al.: Spallation dominated regime 



towards the study of iLezniak & Webbed (Il979l) who showed 
that a truncation of the short path lengths was equivalent to the 
removal of the nearby sources in a DM. As an illustration (not 
displayed here), we find that to obtain the same curve with the 
source hole diffusion model or with the WS model truncated at 
xq - I g cm"^, the size of the hole has to be set around 180 pc 
(when using the median diffusion coefficient). 

As the energy of the events measured in the detectors are 
not well known, the propagation correction cannot be easily 
performed. As explained in the caption of the figure, these 
curves give rough corrections to transform LB -derived source 
abundances to the "more realistic" ones in hole models. For 
illustrative purpose, we apply th i s corre ction to the LB abun- 
dances obtained by iBinns et Figure |9l shows that 
for the specific CR abundances obtained by these authors, the 
LB-propagated abundances (open triangles) display a discrep- 
ancy compared to the SS ones above Z = 60, namely that 
they are twice larger. Previously, a possible enhancement in 
the r-process contribution in this charge range was suggested. 
Applying the approximate correction due to a gas sub-density 
in our diffusion model, a better agreement with SS abundances 
is obtained (filled triangles) for heavier species while fighter 
species are less affected. One must bear in mind that this is a 
very rough estimation. 
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Fig. 9. Circles: Solar System abundances llLoddersI Eool . 

Squares; Measured CR abu ndances combining HEAO and 
Ariel data jBinns et al.lll989l) . The open up triangles represent 
the CR source abundances using a Leaky Box model from the 
above data. Filled down triangles are these values corrected by 
the factor between the LB model to a diffusion model with a 
hole in gas. These factors are taken from Fig.|8l and correspond 
to the median propagation set. For the ranges 44 < Z < 48, 
62 < Z < 69, 72 < Z < 74, the nuclei were considered to be 
pure secondaries. 



5. Conclusion 

The propagation of heavy to UHCRs was considered in a two 
zone diffusion model with various assumptions about the local 
gas density and source distribution. The propagation parame- 
ters used here were previously shown to be compatible with 
B/C and antiproton data, which allows a coherent treatment 
of both heavy and lighter species. Moreover, a previous analo- 
gous study for radioactive species has validated the relevance 
of modelling the observed local sub-density as a circular hole 
of ~ 100 pc in the Galactic disk, increasing confidence in the 
model used here. 

It is well known that for light nuclei. Leaky Box and dif- 
fusion models are equivalent, which is not the case for the ra- 
dioactive species. Using propagation parameters matching B/C 
ratio in a LB and a DM, we explicitly checked that the same 
equivalence exists for the heavy nuclei. However, this equiv- 
alence does not hold when a local sub-density in gas and/or 
sources is taken into account. These new models can be under- 
stood as physical configurations which change the path length 
distribution. They are equivalent to truncation - whose impor- 
tance has previously been underlined for the propagation of 
heavy nuclei - for the case of a hole in sources. To our knowl- 
edge, there is no well-know modified PLD corresponding to the 
case of a hole in gas. 

It was found that for nuclei in the same mass range, these 
models have a weak impact on estimates of CR abundances 
as long as one deals only with pure primaries and pure sec- 
ondaries. The determination of mixed species abundances, e.g. 
Ft, is on the other hand very sensitive to the presence of a local 
gas sub-density. The strength of this effect depends strongly on 
the diffusion coefficient value (and hence on the energy of the 
CRs). 

Whereas this effect is small for propagated nuclei of simi- 
lar masses, it becomes important when considering the whole 
range of nuclei from iron to the actinides. For a typical dif- 
fusion coefficient, it is found that the CR abundances derived 
in a gas hole model are lessened by a factor of 2 compared to 
those evaluated in a LB model. In this case, the discrepancy 
with solar abundances for Z > 60 nuclei is smaller This could 
have important consequences on the of r-process contribution 
needed to explain the measured data. Note that all these results 
were derived independently of the value of the production cross 
sections. 

Further work - using a full propagation network for nu- 
clei - is required to obtain more quantitative results. Before 
this work to be completed, a similar study will first be con- 
ducted for nuclei that are unstable to electronic capture: the 
la tter are known to a ffect source abundances (see e.g.. Table I 
in lLetaw et al1(ll985l) ^ and the local subdensity is likely to have 
an effect (slightly different than the one observed for yS-unstable 
nuclei) on their propagation. 
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Inserting Eq. JA.3l i in the diffusion equation, Eq. iA.H becomes 

d-Ni VcdNi Qi 

' Ni = -2li^6{z) 
K 

(A.4) 



dz^ K dz 



The calculation leading to the solutions for the two zone model 
in cylin drical geometry with out energy losses is presented and 
follows lDoiiato et al.lll2002i) . The diffusion equation describing 
the evolution of the nucleus density N(r, z), including spallation 
and convection due to the galactic wind reads: 

-KAN(r, z) + Vc^+ 2hY6(z)N{r, z) = q{r, z) , , 

oz In the halo, the RHS term of Eq. iAAl is not present. 

where q(r, z) = q{r) x 2h6{z) is the source distribution, K is the Considering the boundary condition NAz = D = 0, the so- 



1 

+ — 
K 



lution in the halo is given by 



diffusion coefficient, Vc is the velocity of the galactic convec- 
tion wind, h is the thickness of the disk, and Y - ncrv is the . I 
spallation reaction rate. The Dirac distribution is needed as the 

spallation and sources are only present in the disk. Considering ' \2K I ' sjnh^^^ 

a local density with a radius a, the density is then given by 

n - nisM®{r - a), where 0(r - a) is the Heavyside function, where 5, is a constant defined as 
The space of the Bessel functions is well adapted to this geom- 
etry and we use the following decompositions, with the i-th 
zero of the Bessel function 7o: 



(A.5) 



■2 \ 



N{r,z)^Y,^i{z)Jo[^i'-^ 

i 

0(r-a)A?(r,z) = 2Q,7o(^,-^) 



(A.l) 

The solution in the thin disk, A^,(0), is obtained by integration 
Eq. iAA\ across the disk between z = -h and z = +/; with 
li 0. Eq. iAA\ becomes 



2Nm - 2^Ni(0) + 2h^- ^2hTBiNi{0) 

K A A 



1 a 

-4hr 



with Qj following 



(A.6) 



J*' 



^i) Jail 



pNip,o)eip--)M(ip)dp 

K 



The continuity between the halo and the disk is established by 
inserting the halo solution (Eq. JA.5» in Eq. iA.bi . Defining A, 

as 



In the space of Bessel functions, the diffusion equation be- 
comes 

Using the property of the Bessel functions 

■ ^[-(iPM^jP)Jii(iP)+ 



Ai = Vc + KSiCothl-^j + 2hYBi 



one finds 



NM = ^2h 
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AYa 
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(A.7) 



(A.8) 



(jpJl{(jp)Jo{(ip)] 

[ V[-/o(^/P) + -/?(^/P)] 



for / 7^ j 
for i = 



reads, 

BiNi for i - j 

with Bi and M,; being respectively 



It appears that to calculate Ni{Q), one needs to know the values 
for all the other orders of the Bessel decomposition. We use a 
perturbative method to compute these quantities. At the zero-r/z 
order, we have A^^'''(0) = 2/1^,/A,, and calculate recursively the 
(« + l)-th order as 



(A.3) A?("+"(0) 



a'' 



4Ya 



, ■ > MijNf\Q) 
qiRjj(^i)j;f ' ^ \ 



(A.9) 



B, = 1 



Mij = 



R^ Aid) 
1 



until the convergence is reached. In practice, the convergence 
is reached quite rapidly, after 5 iterations. Using Eq. \A.\\ . the 
density of nuclei in the physical space is then given by 



„^ is; sinh(5,-^) / r\ 



(A. 10) 
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A.1. Primaries and secondaries 



References 



The previous derivation is valid for any source term q{r). For 
the primaries, we consider the two situations where there are 
sources or not in the hole, but in each case, we assume a con- 
stant source distribution qo with r. In that case, there are ana- 
lytical expressions for the components of q{r) in Bessel space, 
namely 



qir) = 1 



q(r) = 0(r - b) 



1i = 



The secondaries are only produced by spallation of the primary 
nuclei on the ISM, so their source term is q{r) - FpsN'^ir), 
where F,,, is the production reaction of the secondaries and 
N''{r) is the density of primary nuclei. As a consequence, for 
the secondaries we compute Eq. ( IA.8> simply using ^, = F^^A^'^ 
where Nf, referring to the primaries, have been previously de- 
termined. 



Appendix B: The Truncated Weighted Slab 

In this work, we use t he truncated weighted s lab (T WS) ap- 
proach as described in lClinton & Waddingtonl(ll993l) . The es- 
cape mean free path and path length distribution are defined 
respectively as. 



5.5 g cm ^ 

\7.6GV/ 



for R < 7.6 GV 



5.5x1—— forR >7.6GV 



and. 



■ exp 



xo- X 







for x> Xo 
otherwise 



where R is the rigidity and xq the truncation of the shortest path 
lengths (xo=l g cm"- in this work). When neglecting the energy 
losses and assuming only one primary parent p for a secondary 
s, one finds the truncated weighted slab densities of primaries 
and secondaries to be 



-qpe 



(Ap — As)Ap; 



■qp 



{Aesc + Ap) 



(Aesc + As) 



where Ap (resp. As) corresponds to the mean free path of a pri- 
mary (resp. secondary) with respect to its total destruction cross 
section and where Aps is the mean free path of a primary rela- 
tively to the secondary production cross section. 
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